Generalized modified gravity with the second order acceleration equation 
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In the theories of generalized modified gravity, the acceleration equation is generally fourth order. 
So it is hard to analyze the evolution of the Universe. In this paper, we present a class of generalized 
modified gravity theories which have the acceleration equation of second order derivative. Then both 
the cosmic evolution and the weak-field limit of the theories are easily investigated. We find that 
not only the Big-bang singularity problem but also the current cosmic acceleration problem could 
be easily dealt with. 

PACS numbers: 98.80.Cq, 98.65.Dx 



I. MOTIVATION 

In history, the motivation for modifying GR (General 
Relativity) mainly comes from the fact that GR is not 
renormalizable. So it can not be conventionally quan- 
tized. In the first place, Utiyama and DeWitt showed 
that the renormalization at one-loop requires the higher 
order curvature terms in the action of gravity theories 
[l|. Secondly, Stelle showed the corresponding gravity 
theories with these higher order terms are indeed renor- 
malizable Q). Finally, when quantum effects or string 
theory are taken into account, the effective low energy 
gravitational action also requires higher order curvature 
invariants @-|5|. So it was generally believed that the 
modifications to GR would be important only at the 
scales of very close to the Planck lengthy or Planck en- 
ergy. Consequently, both the Big-bang singularity and 
black hole singularity are expected to be absent in the 
modified gravity theories [a-lllj . This is the belief before 
1998. 

However, with the discovery of cosmic acceleration in 
1998 [H, [lj], one realize that GR may also need to be 
modified on very large scale or at very low energy (or 
very weak gravitational field). These constitute the in- 
frared modifications to GR, for example, the GDP (Dvali- 
Gabadadze-Porrati) model [IH, the 1/R modified grav- 
ity model [HI and so on. Here we shall not produce 
an exhaustive list of references on modified gravity, but 
we prefer the readers to read the nice review paper by 
Sotiriou and Faraoni [l6| and the references therein. In 
general, the equations of motion for the generalized mod- 
ified gravity are of fourth order and one can expect the 
particle content of the theory would have eight degrees 
of freedom: two for the massless graviton, one in a scalar 
excitation and five in a ghost-like massive spin two field 
[l7| . The presence of ghost leads one to accept unphysi- 
cal negative energy states in the theory and the property 
of unitarity is lost 18]. This ghost problem is closely 
related to the higher order property of the theories. 

So the purpose of this paper is to seek for the second 



order theories of gravity, at least in the background of 
spatially flat FRW (Friedmann-Robertson- Walker) Uni- 
verse. Except for satisfying the requirement of second- 
order, the theories also meet ghost free conditions. Due 
to the property of second order of acceleration equation, 
the resulting Friedmann equation remains first order and 
the cosmic evolution of the universe is easily deal with. 

The paper is organized as follows. In section II, we 
briefly review the generalized modified gravity theories. 
The equations of motion are presented. In section III, 
we propose the Lagrangian for the generalized modified 
gravity which are both ghost free and second order (in the 
background of spatially fiat FRW Universe). In section 
IV, we investigate the cosmic evolution of some specific 
models of Lagrangian. In section V, we investigate the 
weak field limit of these models. Section IV gives the 
conclusion and discussion. 

We shall use the system of units in which G = c = h = 
1 and the metric signature (— , +, +, +) throughout the 
paper. 



II. GENERALIZED MODIFIED GRAVITY 

The generalized modified gravity theories have the ac- 
tion of the form [19( 

S = J d 4 x^g(^-[R + f(R, P, Q)}+J? m \ (1) 

where / is a general function of the Ricci scalar R and 
two curvature invariants, 



P — R^yR^ , Q = R^axR^ 



(2) 



which are of lowest mass dimension and parity- 
conserving. R^v and R^ va \ are the Ricci tensor and the 
Riemann tensor, respectively. «5f m is the Lagrangian den- 
sity for matters. 
If we define 



f _ df df df 

fR= M , fP=QP> «=flQ' ^ 
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then we obtain the generalized Einstein equations [1J 
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is the generalized Einstein tensor, G M „. The same as 
the Einstein tensor, it satisfies the Bianchi identity 







(5) 



But different from the Einstein tensor which is up to 
second order derivative, 9ffu, is up to fourth-order. So 
Eqs. (TJ| are usually a set of fourth-order differential 
equations except for the case of a cosmological constant 
/ = 2A. The existence of higher derivatives in the 
equation of motions suggests that one would always find 
ghosts in a linearized analysis. Actually, it can be ar- 
gued by considering the Cauchy problem, the gauge sym- 
metries and constraints on the theory (see for instance 
[13]) that the generalized modified gravity will contain at 
most eight degrees of freedom: two for the usual mass- 
less graviton, one for a scalar field and five for ghost-like 
massive spin two excitation. The ghost problem leads 
one to accept negative energy states in the theory. So 
the property of unitarity is lost [l8j |. 

However, Comelli [2(j, Navarro and Acoleyen (2lj 
showed that with a suitable choice of parameters, the 
theory would be ghost-free. Actually, they showed that 
the general Lagrangian of the form Jz? = Jz? (R, P — 4Q) 
are ghost free. So in this case it is left with only an extra 
scalar degree of freedom to the gravitational sector. In 
the next section, we shall seek for the Lagrangian which 
gives the acceleration equation of second order derivative 
and the corresponding theories are ghost free. 



III. THE LAGRANGIAN 

The spatially flat FRW metric is given by 
ds 2 = -dt 2 + a [tf (dr 2 + r 2 dn 2 ) , 



(G) 



where a(t) is the scale factor. Given the metric, we could 
calculate the Ricci scalar R and the curvature invariants, 

P, Q, 



R = 6IH + 2H- 



H + H' 



P = 12 
Q = 12 



H + H' 



H 2 (H 2 + H 



H 4 



H 



(7) 



where H = a/a is the Hubble parameter and dot denotes 
the derivative with respect to the cosmic time t. 

The second derivative of the scale factor, a (in H), 
is present in R, P, Q. If / oc R, the corresponding 
equations of motion are the Einstein equations. They 
are second-order differential equations. In this scenario, 
H appears linearly in the Lagrangian. However, if H 
appears nonlinearly in the Lagrangian just as contributed 
by R 2 , P and Q, the corresponding equations of motion 
would be fourth-order differential equations. 

It is not hard to conjecture that, if we are able to make 
H disappear in the Lagrangian / such that it is uniquely 
the function of Hubble parameter H, the resulting equa- 
tion of motion must be of second-order differential equa- 
tion. Then how to make H disappear? We can examine 
the proper combination of R 2 , P, Q. 

To this end, let us calculate 



aR 2 + fiP + jQ = 12 (3a + /3 + 7) H 2 



So if we let 



namely, 



-12 (12a + 3/3 + 2 7 ) (H 2 H + H 



12a + 3/3 + 2 7 = , 



we would obtain 



7 = -6a - -(3 



(36a + 6/3) H 2 



(8) 
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Is H negative or positive? In order to answer this ques- 
tion, let us resort to the Einstein equation in FRW Uni- 
verse 



H oc — (p + p) 



(12) 



where p, p are the total cosmic energy density and pres- 
sure, respectively. It is apparent H is nonpositive in the 
history of the Universe which mainly covers three epochs 
dominated by radiation, matter and cosmological con- 
stant, respectively. Therefore, we assume 

H < , (13) 

in the following. Thus 

V7= y/-36a- 6/3 (-H\ . (14) 
It is apparent a and j3 should obey 

-36a-6/3>0. (15) 
If we define J as follows 
1 



J 



IR 
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\/ai? 2 + /3P + 7Q 



(16) 
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then we have 



J = H 



in the background of FRW Universe. 

Now we could conclude that, for the general function 
of /(J), the Lagrangian density 



if 



1 

16tt 



[R + f ( J)} + Sf n . 



(17) 



always leads to the second order equation of motion in the 
background of spatially flat FRW Universe. Of course, we 
have assumed the energy-momentum tensor contributed 
by the matters is up to second-order. 

In general, the above theories of gravity would con- 
tain the massive spin two ghost field in addition to the 
usual massless graviton and the massive scalar [l6j]. But 
the .f(R) theories of gravity are found to be ghost free. 
Ref. [lU and Refs. [2C|,|2l| showed that the models given 
by 

Sf=-L[R + f(R,4P-Q)], (18) 
are also ghost free. In view of this point, we should let 

(19) 



£ = -4. 
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Taking account of Eq. (fit)]) , we have 



24 



-a 



(20) 



So J is found to be 



J = 



R + y/6 (4P - Q) - 5i? 2 . (21) 



Using the Gauss-Bonnet invariant: 



G = R-4P + Q 



function of R, P, Q or R,G according to Eq. (|23[) . In the 
background of four dimensional FRW Universe, we have 
J = H. But it is not the case for other spacetimes. 

Here we would like to point out that the above con- 
clusion is valid only for four dimensional FRW Universe. 
If the dimension of spacetime is greater than four, the 
variation of Eq. (|2"4")l would lead to a fourth order ac- 
celeration equation. The reason could be understood as 
follows. In the gauss-Bonnet gravity, if the dimension of 
spacetime is four, the Gauss-Bonnet term turns out to 
be a topological invariant and so it makes no contribu- 
tion to the equation of motion. In higher dimensions, the 
Gauss-Bonnet term is not a topologically invariant and 
so it would contribute to the equation of motion. For the 
modified Gauss-Bonnet gravity, the equation of motion 
is usually fourth order even in the four dimensions. 

Actually, the Gauss-Bonnet term is generalized in the 
Lovelock gravity theory [2J]. It is found that the equa- 
tions of motion in Lovelock gravity are of second or- 
der in any spacetime with any dimensions, needless to 
say in the background of four dimensional FRW Uni- 
verse. However, there is a notable difference between the 
Lovelock gravity and the generalized modified gravity in 
Eq. (|2"4")l . The Lovelock gravity is an ultraviolet modi- 
fication to General relativity and unable to achieve the 
infrared modification to gravity. 



IV. SECOND-ORDER ACCELERATION 
EQUATION AND THE FRIEDMANN EQUATION 

In this section, we shall derive the acceleration equa- 
tion and the Friedmann equation from the Lagrangian 
density Eq. (fT7|) or Eq. (l24l) . In the background of 
spatially flat FRW Universe, the Lagrangian function is 
given by 



L 
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[R+f(J)]+&„ 



(25) 



(22) 



Substituting the Lagrangian function into the Euler- 
Lagrange equation 



we have 



R + \JR? - 6G . 



Then we recognize that the Lagrangian 



l07T 



f(J)} , 



(23) 



(24) 



is actually the modified Gauss-Bonnet gravity which has 
been studied in Ref. [HJ]. In general, the acceleration 
equation in the modified Gauss-Bonnet gravity is four or- 
der. However, for our specific case (Eq. (|2"3"10 . we would 
obtain a second order acceleration equation from the La- 
grangian Eq. (|2"4"]l . We note that, in Eq. (|2"4")l . J is the 



dt 2 



dL 
da 



d_ 

dt 



dL 
da 



dL 
da 



(26) 



we obtain the acceleration equation 



2H + 3H 2 - l -Hf + -/ - l -H.f" = -Sirp , (27) 
2 2 b 

where the prime denotes the derivative with respect to 
H . It is obvious the acceleration equation belongs to the 
second-order differential equations. 

Using the acceleration equation, we could obtain the 
Friedmann equation from the energy-conservation equa- 
tion 
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-£ + 3H(p + p)=0, (28) 

as follows 

3H 2 + -/ - hlf' = 8np . (29) 

The left hand side of the Friedmann equation is only the 
function of Hubble parameter. It is remarkably simple in 
the investigation of the evolution of the Universe. 

Eqs. (f2"TI |2"5I 12"9")) and the equation of state for matters 

p = p( P ) . (30) 

constitute the main equations which govern the evolution 
of the universe. Among the four equations, only three of 
them are independent. For convenience, we always focus 
on Eqs. (EilMllSni)- 

V. SOME EXAMPLES 



When n = 1, we have 

f = l&m 1 H\\\H . (34) 

Now we have the conclusions as follows. 

(1). When / = Wirr/JlnJ (for n = 1), we have the 
following Friedmann equation 

3H 2 = 8tt (p + tjH) . (35) 

It is the same as the Friedmann equation given by the 
DGP modified gravity (26J. The equation can be rewrit- 
ten as 



3H 2 = 87175 + -rj 2 + -rjJrf + 96np . (36) 
6 6 

If we define 



1) = ^Qpi/n/A , (37) 

we have 



In this section, we shall study some specific and inter- 
esting forms of f(H). 



If / = 
Eq. 



ACDM Model 



T671A with A a constant, we obtain from 



3H 2 = 8tt (p + A) 



(31) 



This is the Friedmann equation for ACDM (A-Cold Dark 
Matter) universe. Although the ACDM model provides 
an excellent fit to the wealth of high-precision observa- 
tional data, on the basis of a remarkably small number 
of cosmological parameters [25], it is plagued with the 
well-known cosmological constant problem and the cos- 
mic coincidence problem which prompted cosmologists to 
look for other explanation for the observed accelerated 
expansion. 



B. Power law for / 

Assume the energy density contributed by / is the form 
of a power law r\H n with rj a positive constant and n an 
integer. Then from the Friedmann equation Eq. (|29[) we 
have 



1 



2 J 2 
Thus / is derived as 



/ = 16-7T77 



H n 
n-1 



(32) 



(33) 



3H 1 = 8tt 



EL 
2 



2p- 



2 



(38) 



Here pi is a constant energy density. The pi terms are 
investigated as the candidate of dark energy in many lit- 
eratures, for example, 2jJ and references therein. Dif- 



ferent from the cosmological constant, this dark energy 
density increases with the increasing of background en- 
ergy density p. So the cosmic coincidence problem is 
greatly relaxed. But it is argued that the DGP model is 
disfavored by the history of cosmic structure formation 
(28j because of the fast increasing of dark energy density 
with redshifts. 

(2). When / = — ?yJ~ 2 (for n = —2) with 77 a constant, 
we obtain from Eq. 



3H A 



-r]H~ 2 = 87175 



The above equation can be rewritten as 



3H 2 = 47177 + W 16ttV + -77 . 



If we define 

128 2 2 
V = — 7T Pi , 

then Eq. (|40| can be rewritten as 



3H 2 



Air 



Ap 2 



(39) 



(40) 



(41) 



(42) 



It is apparent pj plays the role of a constant energy den- 
sity. When p ^> pi, it restores to the standard Fried- 
mann equation. When p <C pi, the Universe evolves into 
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a de Sitter phase. We have shown that this model could 
interpret the current acceleration of the Universe [29j . 
Different from the DGP model, the dark energy density 
in this scenario decreases with the increasing of redshifts. 
So the cosmic coincidence problem is also relaxed. 

(3). When / = 16n-q,J 4 /3 (for n = 4), we have the 
Friedmann equation as follows 

3H 2 = 8tt (p + i]H 4 ) . (43) 

It is a quadratic equation of H 2 . Mathematically, we 
would have two roots for H 2 . But physically, only one of 
which could reduce to the standard Friedmann equation 
in the limit of small p. So the physical root takes the 
form of 




0.2 0.4 t 0.6 0.8 1 



3H Z = 



Define 



167T77 



n 



y/9 - 256ir 2 pr) 



128tt 2 Pu 



(44) 



(45) 



with pu some positive constant. The Friedmann equation 
Eq. ((H]) turns out to be 



3H 1 = 8npu 1- ,/l 



2p 
Pu 



(46) 



Here pu plays the role of a constant energy density. It is 
apparent p should obey p < pu- So, to the zeroth order 
of p/pu, we obtain the standard Friedmann equation. To 
the first order of p/ pu, we obtain the Friedmann equation 
in the Randall-Sundrum brane world model 30] 



Putting 



3H 2 = 8tt[p+^- 
2pu 



Pu = 1 , P 



(47) 



(48) 



and using Eq. (|4l)|) , we plot the evolution of the scale fac- 
tor a and the Hubble parameter H in Fig.[TJ respectively. 
It shows that the universe is created in finite time with 
finite scale factor and finite Hubble parameter. So the 
Big-Bang singularity is avoided. 

It is apparent the energy density is also finite from 
Eq. (l4l)|) . Then how about the pressure and its higher- 
derivatives? If they are irregular, some weak singulari- 
ties would appear. We find that there are regular and 
there are no weak singularities. The proof are as follows. 
Taking account of Eq. (|4"51) and the energy conservation 
equation, we obtain the pressure p = p/3. Since p is 
regular, p is also regular. Then with the help of energy 
conservation equation, we find the higher derivatives of 
pressure are regular. 



FIG. 1: The evolution of the scale factor a and the Hubble 
parameter H with respect to the cosmic time t. It shows that 
the universe is created in finite time with finite scale factor 
and finite Hubble parameter. So the Big-Bang singularity is 
avoided. 



C. Past de Sitter universe 

In the last part of subsection B, we find the universe 
could be created in finite time with finite scale factor 
and finite Hubble parameter. So it seems there exists a 
starting point of time. Different from this case, in this 
section, we present a past de Sitter universe. In this 
scenario, the universe starts from a de Sitter phase. So 
the history of the universe is infinite. In other words, 
there does not exist a starting point of cosmic time. 

In order to construct a past de Sitter universe, we could 
explore 



/ = 6J 2 - 4Jy/67rpu tanh -1 



/ 3J 2 
8npu 



(49) 



where pu is a constant energy density. Substituting it 
into Eq. (|29|) , we obtain 



3H 2 



3H 2 , 



m z = 8tt 



1 + p/pu 



(50) 



(51) 



To the zero order of p/pu *C 1, it restores to the Fried- 
mann equation in General Relativity. To the first order 
of p/ pu, we have 



3H 2 



8tt[p-S- 
Pu 



(52) 



It is the Friedmann equation in the theory of loop quan- 
tum gravity [3l|, HH 



6 



On the other hand, when p/ pu S> 1, we obtain 



3H 2 = 8irpi 



(53) 



This is for de Sitter universe. Taking into account all 
matter sources which include relativistic matter (radia- 
tion), baryon matter, dark matter and the cosmological 
constant, we can rewrite the Friedmann equation Eq. (|51[) 
as follows 



f2 r0 



(54) 



Observations show that f2 r o = 8.1 • 10 -5 , f2 m o = 0.27, 
fl\ = 0.73, which are the ratio of energy density for radi- 
ation, matter (including baryon matter and dark matter) 
and cosmological constant in the present-day universe. 
The dimensionless Hubble parameter h is defined by 



(55) 



with Hq the present-day Hubble parameter. Qu is the 
dimensionless energy density contributed by pu- 

Taking VLjj = 10 123 (which represents pu is Planck 
energy density), we plot the evolution of the re-scaled 
dimensionless Hubble parameter M = log 10 h with re- 
spect to the re-scaled scale factor N = In a in Fig. @. 
There are three epochs A, B, C in the total life of the 
Universe. The epoch of A corresponds to the de Sitter 
phase. The Universe exponentially expands (inflating) in 
this period. Then the inflation stops around the redshift 
of z ~ 10 30 . The epoch of B is dominated firstly by the 
radiation and then by the matter. It stops around the 
redshift of z ~ 0. The epoch of C is dominated by a small 
cosmological constant. It is the future de Sitter phase. 



VI. WEAK FIELD LIMIT 

In this section, we shall study the weak field limit of the 
above modified gravity theories. Following Ref. [20l [2ljj . 
we expand the action in powers of the curvature pertur- 
bations. Then it can be shown that at the bilinear level 
the linearisation of the theory over a maximally symmet- 
ric spacetime will be the same as the theory obtained in 
17, 331. 




2ml 



where C^vaX is the Weyl tensor and we have defined 




-100 -80 -60 -40 -20 
N 



20 40 



FIG. 2: There are three epochs A, B, C in the total life of the 
Universe. The epoch of A corresponds to the de Sitter phase. 
The Universe exponentially expands (inflating) in this period. 
Then the inflation stops around the redshift of z ~ 10 30 . The 
epoch of B is dominated firstly by the radiation and then by 
the matter. It stops around the redshift of z ~ 0. The epoch 
of C is dominated by a small cosmological constant. C is the 
future de Sitter phase. 



A = </ - Rf R + R 2 Urr/2 - fp/4 - f Q /6) 

+R 3 (/hp/2 + f RQ /3) 

+i? 4 (/pp/8 + /Q Q /18 + /p Q /6)> , (57) 
6 = <fn-RfRR-R 2 (fR.p + 2f RQ /3) 

-R 3 (/pp/4 + /qq/9 + /pq/3) > , (58) 
- 2 = <3f RR + 2fp + 2f Q + R(3f RP + 2f RQ ) 

+R 2 (3/pp/4 + f QQ /3 + f PQ ) > , (59) 



-2 _ 



= -</p+4/q >o 



(60) 



Here < 



>o represent the values of the correspond- 



ing quantities in the background of some spacetime and 
Frr = e ^ c - It is apparent for the action, Eq. (|24l) . 
the inverse of mass squared of the ghost is m^ 2 = 0. 
Thus there is no ghost in the spectrum. But there is still 
an extra scalar with the mass rriQ. The A term behaves 
as the vacuum energy and the 8 term contributes the 
variation of the gravitational constant, 



SG 
~G 



1 



(61) 



In TABLE. I and TABLE. II, we calculate the parame- 
ters for the five models: 



(a) 
(b) 
(c) 
(d) 

(e) 



/ = -16ttA; 

/ = 167T77</ln J ; 

f = -vJ- 2 ; 

/ = 16tt77 J 4 /3 ; 

/ = 6 J 2 - 4J A /B7rp^tanh" 



V ^Pu ) ' 
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Models 


A 


5 


777,,, ^ 


777,^ Z 


(a) 


A 

A 


U 


U 


u 


(b) 


-^/A~p7(9 In 2 + 3 In 7T + 3 In A - 8 - 3 In 3)/16 


- ^/p//A(-3 In 3 + 4 + 9 In 2 + 3 In vr + 3 In A)/32 


mo = 





(c) 


P?/A 


-tf/(4A a ) 


mo = 





(d) 








m = 





(e) 








m = 






TABLE I: The parameters of five models in the background of de Sitter spacetime. 



Models 


A 


6 


m " 




(a) 














(b) 








mo = 





(c) 


+oo 


— 00 


mo = 





(d) 








mo = 





(e) 








mo = 






TABLE II: The parameters of five models in the background of Minkowski spacetime. 



in the background of de Sitter spacetime and Minkowski 
spacetime (by taking the limit of A — > 0). From TA- 
BLE I, we see the models a, b, c contribute the non- 
vanishing cosmological constant terms. For models d and 
e, the A terms are zero because they are the ultraviolet 
modification to GR. In the column of 6, we see the models 
6, c contribute the non- vanishing gravitational constant 
terms. Except for the model of a, there exists a scalar 
degree of freedom with vanishing mass in b, c, d, e. 

In the background of Minkowski spacetime, we see 
from the TABLE. II that, except for the model of c, 
the other models make no contribution to vacuum en- 
ergy and gravitational constant. The reason for this is 
that the model c is essentially a modification with inverse 
curvature invariants. So the Minkowski spacetime does 
not solve the corresponding equations of motion. 

VII. CONCLUSION AND DISCUSSION 

In the theories of generalized modified gravity, the ac- 
celeration equation is generally fourth order. So it is 



hard to analyze the evolution of the Universe |15[. On 
the other hand, these theories are also plagued with the 
ghost problem. So the property of unitary of the theory 
is lost [18j. In view of this point, we present a class of 
generalized modified gravity theories which have the ac- 
celeration equation of second order derivative and ghost 
free. Then we explore some specific examples for the La- 
grangian function. We find both the cosmic evolution 
and the weak-field limit of the theories are easily inves- 
tigated. Furthermore, not only the Big-bang singularity 
problem but also the current cosmic acceleration problem 
could be easily dealt with. 
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